Mon. Not. R. Astron. Soc. OOP. [TlfTOl (20121 Printed 3 March 2013 (MN 1^1^ style file v2.2) 



Is there a Size Difference between Red and Blue Globular 
Clusters? 
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ABSTRACT 

Blue (metal-poor) globular clusters are observed to have half-light radii that are 
~20% larger than their red (metal-rich) counterparts. The origin of this enhancement 
is not clear and differences in either the luminosity function or in the actual size of 
the clusters have been proposed. I analyze a set of dynamically self-consistent Monte 
Carlo globular cluster simulations to determine the origin of this enhancement. I find 
that my simulated blue clusters have larger half-light radii due to differences in the 
luminosity functions of metal-poor and metal-rich stars. I find that the blue clusters 
can also be physically larger, but only if they have a substantial number of black holes 
heating their central regions. In this case the difference between half-light radii is 
significantly larger than observed. I conclude that the observed difference in half-light 
radii between red and blue globular clusters is due to differences in their luminosity 
functions and that half-light radius is not a reliable proxy for cluster size. 

Key words: galaxies: star clusters: general - globular clusters: general - stellar 
dynamics - methods: N-body 



1 INTRODUCTION 

Globular cluster (GC) systems are common in disk and el- 
liptical galaxies and in the past decade significant obser- 
vational progress has been made in understanding them. 
These globular clusters are frequently used to help under- 
stand the formation and evolution of their host galaxies. 
However, the observational properties of globular clusters 
are not fully understood. A case in point is the half-light ra- 
dius (r/ii). Several models show that the the half-light radius 
of a cluster should remain fairly constant during its lifetime 
(jSpitzer fc Thuanlll972l : lAarseth fc HeggieHlQQsf l so it is of- 
ten used to compare the structures of globular clusters of 
different ages. Several observational studies, however, show 
that metal-poor (blue) globular clusters have half-light radii 
that are systematically larger (by ~ 20%) than their metal- 



rich (red) counterparts 


(iKundu & Whitmore 19981. 2001 ; 


Kundu et al.l iggg*: 'Puzia ct al." 'l999; 'Larson ct al.' 2001; 


Larsen. Forbes & Brodic 


,2001.; .Barmbv. Holland & Huchra 


2OO2I: Harris et al. 20021; 


Jordan et al. 20041; HarrisI 20091). 



It is not clear what the origin of this discrepancy is 
nor is it certain that it truly represents a difference in the 
sizes o f the clusters as measu red by their half-mass radii 
{rum). iLarsen fc Brodiel (|2003l ) proposed that the observed 
enhancement in the size of the blue clusters could be due to a 
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projection effect. They note both that blue and red globular 
clusters follow different radial distributions and that in the 
Milky Way there is a rel ationship between cluster size and 
galac tocentric distance (|van den Bergh. Morbey fc Pazden 
Il99lf ). They argue that if such a relationship exists in 
all galaxies, the differing radial distributions of red and 
blue clusters could accoun t for the obs erved difference in 
sizes. A detailed survey by iHarrisI (|2009l ). however, showed 
that the ratio of rhi between blue and red clusters does 
not depend on galactocentric distan ce, as would be pre- 
dicted by t he Larson fc Brodiel (|20Q3l ) model. Furthermore, 
iLarsen fc Brodic (,200l ) give no reason why a galactocen- 
tric distance-globular cluster size trend should exist in other 
galaxies. 

IJordanI (|2004l ) explained the difference as a result of dif- 
fering stellar evolution processes in metal-poor and metal- 
rich stellar populations. Metal poor stars lose less mass and 
have longer main-sequence lifetimes than metal-rich stars. 
Assuming that the difference between red and blue stellar 
populations is metallicity, this leads to differing luminosity 
functions in blue and red GCs. Using multi-mass Michie- 
King models with fixed half-ma ss radii and a stellar popu- 
lation with an age of 13 Gyrs, IJordanI (|2004l ) was able to 
re-produce the observed rni enhancement in blue globular 
clusters. However, these models are not dynamical simula- 
tions of globular cluster evolution and the results are valid 
only if rhm is the same in red and blue globular clusters. 
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By contrast [Schulman. Glebbeek fc s53 (|2012l ) used di- 
rect A'^-body simulations to investigate the differences in dy- 
namical evolution between blue and red open clusters. They 
also assumed that the colour of a cluster reflects its metal- 
licity and argued that, because metal-poor stars lose mass 
more slowly than metal-rich stars, blue clusters will lose less 
mass to stellar evolution over their lifetimes than will red 
clusters. This reduces the gravitational potential of red clus- 
ters, causing them to expand and become larger than blue 
clusters. However, scattering interactions between the more 
massive metal-poor stars will be more energetic than those 
between less massive metal-rich stars and the blue clusters 
will experience stronger two-body heating. This will cause 
them to expand relative to the red clusters once the initial 
phase of rapid mass-loss is concluded, ff the effect of two- 
body heating is strong er, blue clusters will be physically 
larger than red clusters. ISchulman. Glebbeek fc Sillj (|2012l ) 
showed that the half-mass radii of their simulated blue clus- 
ters were indeed ~ 20% larger after several initial half-mass 
relaxation times {trh) than the half-mass radii of their simu- 
lated red clusters. They also found little difference between 
the ratio of rhm and rni between blue and red clusters. 
Therefore they conclude that the enhancement in rhi in blue 
globular clusters is due to an actual size enhancement and 
has a dynamical origin. The ISchulman. Glebbeek fc Sillj 
models, however, contain 10-100 times fewer stars 
than are present in globular clusters and consequently re- 
lax and dissolve when they are only a few hundred Myr 
old, far younger than the 10-13 Gyr age of many globular 
cluster systems. The relationship between stellar evolution 
timescales and dynamical timescales is also quite different 
in these simulations than in GCs. 

In this paper I re-visit the problem of size differences 
between blue and red globular clusters by analyzing a set of 
Monte Carlo star cluster models. These are self-consistent 
dynamical simulations that contain a similar number of stars 
to globular clusters, include parametrized stellar evolution 
and primordial binaries. I will investigate whether my blue 
globular clusters are larger than the r ed on es, as was re- 
ported by ISchulman. Glebbeek fc Silld \2Qli ). If they are, I 
will determine if rni is enhanced to the s ame d egree as rhm 
and if the processes reported bv lJordAnI (|2004l ) have an sig- 
nificant effect. I will also be able to determine weather or 
not Thi is a good observational proxy for rhm- 



2 NUMERICAL MODELS 

The simulations in this paper were performed using a Monte 
Carlo method to self-consistently simulate the dynamical 
evolution of a star cluster in the Fokker-Planck two-body 
relaxation limit. The primary advantage to this method 
over direct A'^-body is speed; the Monte Carlo code scales 
with 0{N^) - 0{N'^), where is the number of stars in 
the system, as opposed to 0{N^) - 0{N'^) for direct iV- 
body codes. Thus it it is possible to run globular cluster- 
sized simulations over a full Hubble time of dynamical 
evolution, a task currently impossible to do using a di- 
rect A''-body code. Unlike standard Fokker-Planck codes, 
however, the Monte Carlo method treats the cluster as 
an ensemble of particles and can provide the same star- 
by-star information as a direct Ai'-body simulation. This 



Table 1. Parameters of the simulations. The first column gives 
the identifying label, the second gives the metallicity, the third the 
total mass of the cluster, the fourth the initial half-mass radius 
and the fifth the initial half-mass relaxation time. 

Simulation Parameters 



Simulation 


Z 


M [Mq] 


rhm [pc] 


trh [Myr] 


red21 


0.02 


3.61 X 10^ 


7.14 


3.54 X 10^ 


red37 


0.02 


3.63 X 10^ 


4.05 


1.51 X 10^ 


red75 


0.02 


3.62 X 10^ 


2.00 


5.25 X 10^ 


blue21 


0.001 


3.60 X 10^ 


7.14 


3.55 X 10=* 


blue37 


0.001 


3.62 X 10^ 


4.05 


1.51 X 10^ 


blue75 


0.001 


3.62 X 10^ 


2.00 


5.25 X 10^ 



also makes it much simple to include prescriptions for stel- 
lar evolution and strong few-body inte raction s. The Monte 
Carlo code is described in detail by ICiersj (|l998l . l200ll . 
l2006h while the stro ng few-body interactions are calcu- 
lated as described in iGiersz fc SpurzemI (|2003l l. The code 
computes parametriz ed evolution for each sta r and bi- 
jHurlev. Pols fc Toutllioool ) and BSE 



nary using the SSE 



jHurlev. Tout fc Polsl 



2002h ste llar evolution prescriptions 



iGiersz. Heggie fc HurlevI (|2008l ). These are essentially the 



same prescriptions as used bv Ischulman. Glebbeek fc Silld 
(|2012h . The code has been compared to, and provides excel- 
lent agree ment with, both direct N-body simulations and ob- 
serva 1 :ions llCiersz. Heggie fc HurlevI l2008l: iHeggie fc Giersd 



l2008l : IGiersz fc Heggie(l2009l : iHeggie fc Giersj|2009i r 

The simulations used in this paper are a subset of 
a collection of simulations that were originally performed 
to investigate the dynamical creation of black hole bina- 
ries i n globular clust e rs. T hey are described in some de- 
tail in iDowning et~al] (|2010l '). All simulations are initialized 
as 5 X 10^ particle Plummer spheres with 10% primordial 
binaries. Initial binary parameters are chosen acco rding to 
the ei genvalue evolution and feeding algorithms of iKroupal 
(|l995h . The initial stellar density is controlled by the ratio 
of the tidal radius (rt) to rhm- For all simulations rt = 150 
pc and I use rt/rhm G {21,37,75}, corresponding to initial 
stellar number densities within rhm of ~ lO'^ , 10"* and 10* re- 
spectively. Th e masses are between 0.1 and 150 Mq and are 
drawn from a iKroupa. Tout fc Gilmor3 (|l993D initial mass 
function with a low-mass slope of q; = 1.3, a high mass 
slope of Uh ~ 2.3 and a break mass of Mbreak ~ 0.5 Mq. 

The simulations come in two metallicities, Z = 0.02 
and Z = 0.001. Throughout the paper I will refer to the 
Z = 0.02 clusters as the "red" clusters and the Z = 0.001 
clusters as the "blue" clusters. One significant difference be- 
tween the red and blue simulations is the treatment of black 
holes (BHs). Metal-poor stars lose less mass to line-driven 
winds than do metal-rich st ars and ex perience more mat- 
ter fallback after supernovae jBelczvnsk i. Kalo gera fc Buli^ 
l2002l : iBelczvnski et aLlbOOel ). This produces both more and 
more massive BHs in the blue clusters. BHs also receive 
natal kicks drawn from a Maxwellian distri bution with a 
peak at 190 km/s (|Hansen fc PhinnevI 119971 ) . significantly 
higher than the escape velocity of a GC. This kick is then 
reduced according to the amount of fallback on to the BH 
during the supernova (|Belczvnski. Kalogera fc Bulikll2002l ). 
Because the amount of fallback is larger in the blue clusters, 
the BHs in these clusters receive a larger kick reduction and 
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Figure 1. Ratio of ri^i and rji^ in blue globular clusters as com- 
pared to red globular clusters for clusters with different initial 
concentrations. The top plot gives the ratio after 13 Gyr, the 
middle the ratio after three initial half-mass relaxation times, and 
the bottom the ratio after one initial half-mass relaxation time. 
One sigma error bars are given. Half-mass and half-light radii are 
clearly larger in blue clusters. 



are less likely to escape. The combination of the larger BH 
masses and higher BH retention rates leads to significant 
differences in the BH populations of blue and red clusters. 

The combination of initial concentration and metallicity 
give 6 different sets of initial conditions, each of which is then 
independently realized 10 times to constrain the stochastic- 
ity inherent in collisional stellar dynamics. All values in this 
paper are an average over all ten realizations unless specif- 
ically noted otherwise. This yields a set of 60 simulations. 
The basic parameters are given in Table [1] The init ial half- 
mass relaxation times are calculated according to ISpitzej 
(|l987h : 
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Figure 2. The evolution of Rhm,br ^ a- function of time for all 
three different initial concentrations. 



where (m) is the average mass of stars in the cluster. I note 
that each of these simulations takes only ~ 1 day to run 
as opposed to the months needed for direct Ai'-body simula- 
tions. It is this advantage in speed that makes this investiga- 
tion possible. I have performed several further simulations to 
investigate specific effects and these will be described later. 



3 THE RELATIVE SIZES OF RED AND BLUE 
GCS 

In figure[T]I show the ratios of rhm and rhi between blue and 
red GCs. r^i is calculated by projecting the luminosity of all 
stars into annular rings and summing the total luminosity in 
each ring radially outwards until finding the radius contain- 
ing half the total luminosity. The ratio of the value of rhm 
in blue over the value of rhm in red clusters will hear after 
be referred to as Rhm,br while the ratio of the value of rhi in 
blue over the value of rhi in red cluster will be referred to as 
Rhi,br- It is clear that both rhm and rhi are enhanced in blue 
clusters. After 13 Gyr the value of rhm is ~ 20% greater in 
blue clusters, in ro ugh agreement with the value foun d af- 
ter ~ 250 Myr by ISchulman. Glebbeek fc SiUd (|2012l ). By 
contrast the value of rhi is some ~ 40-50% greater, a larger 
enhancement than observed in real clusters. The difference 
in Rhm,br and Rhi,br indicates that the enhancement of rhi in 
blue GCs is not simply the result of a difference in size. That 
Rhi,br is larger in my simulations than in the observations 
suggests that my simulations over-predict the difference be- 
tween blue and red GCs. 

The time evolution of Rhm,br, shown in fig- 
ure [2] is rather similar to th e behaviour found by 
ISchulman. Glebbeek fc Sills! (|2012l ') - rhm starts off smaller 
in the blue clusters but then grows rapidly so the blue clus- 
ters becom e larger than the red ones within an initial relax- 
ation time. ISchulman. Glebbeek fc Sills! (|2012l ) suggest that 
this pattern refiects stronger two-body relaxation between 
the more massive stars in the blue clusters. Both clusters 
experience strong initial mass-loss due to the rapid evolu- 
tion of massive stars and because this effect is stronger in 
the red clusters they initially expand faster. However, the 
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Figure 3. The evolution of the BH population in blue and red 
GCs with an initial concentration of rt/r^^ = 37. The top panel 
gives the number of single BHs, the middle panel the number of 
BH-BH binaries and the bottom panel the average mass of single 
BHs. 



more energetic two-body interactions between the massive 
stars in the blue clusters generate more dynamical heat and 
cause them to expand faster than the red clusters once the 
initial phase of rapid mass-loss is complete. Thus after the 
first few hundred Myrs the blue clusters become larger than 

the red ones. 

ISchuIman. Glebbeek fc Sild (|2012l ) do not determine 
whether the dominant source of the two-body heating in 
the blue clusters is the interactions between a large number 
of stars, all of which are only slightly more massive than the 
stars in the red clusters, or if a small population of very mas- 
sive stars is responsible for the effect. A comparison of the 
time evolution of Rhm,br in figure [2] for different initial con- 
centrations provides a clue. For the lowest density (dynam- 
ically youngest) clusters the value of Rhm,br slowly grows 
to ~ 1.2 and then remains constant whereas in the dens- 
est (dynamically oldest) cluster Rhm,br grows rapidly but 
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Figure 4. The BH and BH-BH binary population in clusters with 
^'t/rhm = 75. The top panel gives the number of single BHs while 
the bottom panel gives the number of BH-BH binaries. 
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Figure 5. A comparison of ratios of and between blue 
and red clusters with rt/ri^^ = 75. 



then begins to decline. This indicates that the densest blue 
clusters start to re-contract compared to the red clusters be- 
cause the population of stars responsible for the expansion 
has been lost. Such behaviour suggests a small population 
of massive objects that can pump a large amount of kinetic 
energy into the blue clusters through interactions but are 
also scattered out of the cluster and depleted. 

BHs are excellent candidates for this population: after 
the first few 100 Myrs they are significantly more massive 
than the average mass in the cluster; they are expected to 
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Figure 6. The same as for figurc[2lbut for clusters whore all BHs 
have been removed by giving them 1000 km/s natal kicks. 




mass - segreRate to the core, become S pitzer-unstable (|Spitzeij 
ll987l : ISigurd sson fc Hernauistlll993l ). and interact strongly; 
they are also fe w in number so they c an be rapidly depleted 
by interactions. iMacKev etlj] (|2008l ') presented simulations 
that show a population of BHs can lead to a significant 
growth in core radius compared to simulations without BHs. 
Figure 8 of the sam e paper suggest s the s ame may be true 
for rhm- Specifically. iMacKev et al] (|2008l ) propose that the 
BHs mass-segregate to the centres of the clus ters where they 
strong ly interact and form BH-BH binaries. iDowning et al.l 
l|201(l ) show that this also occurs in my simulations. These 
binaries scatter single BHs from the cluster centre to its 
outer regions. From there these BHs sink back to the centre 
of the cluster due to dynamical friction. This process does 
'mBH\<i>\ work on the cluster per BH where rriBH is the mass 
of the BH and <j) is the gravitational potential of the clus- 
ter. This work dynamically heats the cluster and causes it 
to expand. 

As discussed in § [21 there can be major differences be- 
tween the BH populations in my red and blue simulations. I 
show the population of BHs in the rt /rhm ~ 37 simulations 
in figure O I choose to analyze this set of simulations be- 
cause there are the least extreme set of initial conditions and 
should be the most generally representative. Figure [3] con- 
firms that the blue clusters indeed have a larger population 
of BHs and that these BHs are significantly more massive 
than those in the red clusters. The blue clusters also produce 
slightly larger numbers of BH-BH binaries slightly earlier 
than the red clusters but the di fference is not nearly a s great. 
This is in line with the model of lMacKev" et all (|2008l 'l where 
the BH-BH binaries do not interact directly with the rest of 
the cluster but only with other BHs. It is these scattered 
single BHs that are responsible for the dynamical heating of 
the luminous stellar population. Although the red clusters 
have not insubstantial BH-BH binary populations compared 
to the blue clusters, they do not have a sufficient number of 
sufficiently massive single BHs to heat the rest of the cluster 
and cause it to expand. 

The BH hypothesis also explains the evolution of the 
most concentrated clusters. Figure [4] shows that the both 
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Figure 7. Mass-loss and number of escapers in blue vs. red clus- 
ters with BHs and with rt/rh^ = 37. The top panel gives the 
cumulative mass-loss due to stellar evolution and escapers, the 
middle the cumulative number of escapers and the bottom the 
total mass as a function of time. 



BHs and BH-BH binaries are rapidly depleted in the 
rt/rhm = 75 clusters. Consequently these clusters loses this 
source of dynamical heat. The Monte Carlo models remain 
in virial equilibrium to within a fraction of a per cent by con- 
struction so this loss of heat must be compensated for by a 
change in the structure of the cluster. The virial radius. 



Rv — 



2K 



K = -Ma'^ 
2 



(2) 



where M is the total mass of the cluster, a the velocity 
dispersion and K the kinetic energy, is a good measure of 
this change. Figure [S] shows that the ratio of i?„ between the 
blue and red clusters increases while there are many BHs and 
BH-BH binaries and then decreases when these objects are 
depleted. Indeed the drop in the ratio of Rv and a sharp 
drop in the number of BH-BH binaries both occur at ~ 6 
Gyr. The ratio of _R„ also tracks the evolution of Rhm,br very 
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Figure 8. Position of last interaction before escape relative to the cluster centre for escapcrs from a single red (left) and single blue 
(right) cluster with rt/rj^^ = 37 as a function of time. The black line is the 10% Lagrangian radius. 



closely. This strongly supports the case that heating by BHs 
is responsible for the size differences between these clusters. 

To confirm that the difference in rhm between my blue 
and red clusters is a result of the different BH populations 
I have run a set of simulations without BHs. For each set 
of initial conditions listed in table [T] I have performed an 
additional five simulations, identical up to and including 
the random seed, but where all BHs are given a kick of 
1000 km/s upon formation. This instantaneously removes 
all black holes while allowing the cluster to continue its evo- 
lution. Values of Rhm.br for the clusters without BHs are 
given in figure [S] and can be compared to the results for 
clusters with BHs given in figure [S] The results are dra- 
matic - the enhancement in rmh of blue clusters completely 
disappears! Indeed the trend observed in figure[2]is reversed: 
the blue clusters are slightly larger at early times - probably 
due to the fact that BHs are more numerous and massive 
in these clusters, resulting in more mass loss when they are 
ejected - and slightly smaller at late times. The effect in 
figure [6] is only ~ 2 — 4%, much smaller than the 20% 
differences observed in figure [21 Because the only difference 
between these two sets of simulations is the BH population, 
this proves that the enhancement in Vhm observed in my 
first set of blue clusters is due to the dynamical activity of 
BHs. 

One effect I observe in my s imulations wit h BHs that 
was not extensively discussed bv iMacKev et all (|2008l ) IS a 
difference in the total escap e rate of stars from the blue and 
red clusters (|MacKev et al.l (|2008l ) examine only the escape 
rate of BHs). In figure [7] I show that the blue clusters lose 
more stars, both by number and by mass, to escape than 
do the red clusters. By contrast the red clusters lose slightly 
more mass to stellar evolution than do the blue clusters. 
The difference in mass lost to escapers is larger and the 
blue clusters lose more mass overall than do the red clusters. 



Figure[S]shows that the majority of these additional escapers 
come from the central regions of the cluster, which suggests 
they are stars that have gained additional energy through 
interactions with the BH sub-system. This additional mass- 
loss from the cluster centre will cause the inner Lagrangian 
radii to expand compared to the overall radius of the cluster 
and enhance the expansion of rhm- Figure [S] gives the total 
energy of each escaping star — Ki, -\- mi,cj> where Ki, 
is the kinetic energy of the star and m* is its mass) as a 
function of the position of its last interaction before escaping 
from the cluster. Most of the additional escapers in the inner 
region do not have particularly high energies. Therefore they 
are unlikely to have been ejected in an interaction with a 
binary. Rather they have gained their energy due to strong 
two-body relaxation in the inner regio ns of the cluste r . This 
is fully consistent with the results of iMacKev et al.l (|2008l ) 
who found that the BH-BH binaries do not interact with 
other cluster stars directly but only with the single BHs. It 
is two-body heating from these single BHs that produce the 
expansion of the cluster and the higher escape rate. 

In figure [TOl I compare the mass-loss and escaper rates 
in blue and red clusters without BHs. The difference in total 
number of escapers almost completely disappears. There is 
still a difference between the mass lost due to escapers but 
the majority of this occurs early and can be attributed to 
the instantaneous ejection of BHs from the clusters. Because 
the blue clusters form more BHs and these BHs are more 
massive, it follows that the blue clusters will lose more mass 
when they are ejected. This may also explain the slightly 
larger values of rhm in blue clusters without BHs at early 
times. This difference has little effect on the total mass of 
the clusters and, unlike the case with BHs, the blue and red 
clusters have approximately the same mass at 13 Gyrs of 
age. 
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Figure 9. The total energy of escaping stars as a function of 
the position of their last interaction before escape for a single 
red (top) and single blue (bottom) cluster with rtlvy^^ = 37. 
The total energy of each star is given in units of the core kinetic 
energy of the cluster. The low-energy peak is primarily ejections 
due to two-body relaxation while the high-energy peak is ejections 
due to interactions with binaries. 
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4 THE HALF-LIGHT RADII OF RED AND 
BLUE GCS 

Figure [T] shows that the value of Rhi,br is larger both than 
Rhm,br and than the value that is observed in real GCs. 
Therefore it is clear from these simulations that Rhi,br is 
not necessarily a good predictor of Rhm.br- To understand 
the difference between Rhi,br and Rhm.br I compare the time 
evolution of Rhi.br in clusters with and without BH in fig- 
ure [TT] The half-light radii are clearly larger in blue clusters 
for all sets of initial conditions, regardless of whether Vhm is 
larger or not. The enhancement in rhi in the clusters with- 
out BHs is smaller but is actually in much better agreement 
with the observed enhancement of ~ 20%, at least for the 
two most initially concentrated clusters. The difference be- 
tween Vhm in the blue and red clusters with BHs contributes 
to the difference in rhi but clearly there are other processes 
at work. 

In figure [121 1 compare the fraction of the total luminos- 
ity in the most massive 10% of stars for the blue and red 



Figure 10. The same as for figurclTlbut for clusters without BHs. 



clusters while in figure [13] I compare the mass-segregation 
in these clusters using the average masses within selected 
Lagrangian radii. The red clusters have a steeper luminos- 
ity function with a significantly larger proportion of their 
luminosity concentrated in their most massive stars. The 
average stellar mass is higher for the inner Lagrangian radii 
in both blue and red clusters so both sets of clusters are 
mass-segregated. For a given Lagrangian radius, however, 
the average mass is similar in both the blue and the red 
clusters. Taking figures [121 and figure[T3]together, this means 
that the red clusters will have a larger fraction of their lumi- 
nosity located within their innermost Lagrangian radii than 
will the blue clusters and thus they will have a smaller value 
of rhi compared to rhm- Consequently the value of Rhi,br 
will be larger than the value of Rhm,br- This explains why 
Rhi,br > 1 in clusters of the same size according to rhm- The 
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Figure 11. The ratio of r^i between blue and red clusters for all 
initial concentrations with and without BHs. The top panel gives 
the values for clusters with BHs while the bottom panel gives the 
values for clusters without BHs. 
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Figure 12. The fraction of the total luminosity contained within 
the most massive 10% of stars for blue and red clusters without 
BHs and with n/r^^ = 37. 



effect will be further enhanced if the blue clusters are also 
larger than the red clusters and explains the greater value 
of Rhi,br in clusters with BHs. However it is not necessary 
for a cluster with a larger value of Vhi to have a larger value 
of riim- This suggests an additional corollary: rhi cannot be 
used to measure the relative sizes of GCs, at least not in a 
straightforward way. 
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Figure 13. From top to bottom: the average mass within the 
10%, 50% and 100% Lagrangian radii as a function of time for 
red and blue clusters without BHs and with rt /r^^ = 37. 



5 DISCUSSION 

My simulations show that blue GCs can be larger than 
red GCs but only if they have a substantial population of 
BHs. Blue and red clusters without BHs will be roughly the 
same size but the red clusters will still have smaller val- 
ues of Vhi due to their steeper luminosity functions. I find 
that the enhancement in rni for blue clusters with BHs is 
significantly larger than that observed value. In the clus- 
ters without BHs the enhancement in rni roughly matches 
the observations, at least for the clusters with the highest 
initial concen trations. This finding supports the results of 
Ijordani (|2004l 'l - that the difference in half-light radii be- 
tween blue and red GCs is the result of different luminosity 
functions in clusters with similar size s - ra ther than the 
resuhs of ISchulman. Glebbeek fc Sillsl (|2012l ). It may also 
indicate that there is little difference between the BH popu- 
lation in blue and red GCs. This will have consequences for 
understanding effect of metallicity on the late-time evolution 

and supernovae of mass ive stars. 

The simulations of ISchulmam Glebbeek fc Sillsl (|2012D 
used the initial mass function of Kroupal ~( 2001 ) , which has 
a high-mass slope of ~ —2.3, and used a maximum Mass 
of 50 M0 (Glebbeek, private communication). Such a mass 
function will yield few BHs so a question remains as to 
why my simulations w ith BHs behave in a qualitatively 
similar way to those of ISchulman. Glebbeek fc Sii5 (|20ia ) 
while the ones without BHs do not. I speculate that the 
similarity is a result of the differing relationship between 
the stellar evolution and the relaxation timescales. In my 
simulations trh is significantly longer than the evolution 
timescale for massive stars, even in the most concentrated 
clusters. Thus by the time relaxation becomes important 
for the structure of the cluster the only remaining popu- 
lation of stars significantly more massive than average are 
the BHs. B y contrast the relaxation timescale s of the sim- 
ulations of Ischulman. Glebbeek fc Silld (|2012h are ~10-50 
Myr, shorter than the evolution timescales of moderately 
massive stars. Therefore these stars remain when two-body 
relaxation starts to drive cluster evolution. Furthermore, the 
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Figure 14. The ratio of r^i to rhm for all simulated GCs. Tiiin 
lines are red GCs while thick lines are blue GCs. The top panel 
gives the values for GCs with BHs while the bottom panel gives 
values for GCs without BHs. 



short relaxation time will allow these stars to mass-segregate 
before they lose a significant amount of their mass to winds 
and nuclear burning. Interactions between these stars will be 
able to act as a heat source in a similar way to the BH sub- 
system in my simulations. I hypothesize that small clusters, 
with relaxation times shorter than the evolution timescale 
for massive stars, have a variety of objects massive enough 
to drive cluster expansion through two-body heating while 
in GCs, where the relaxation timescale is much longer than 
the stellar evolution timescale, BHs are the only remaining 
objects massive enough to produce this effect. 

It is also apparent that the value of Rhi,br is not the 
same as that of Rhm,br and thus it cannot be used as a direct 
estimate of the relative sizes of GCs. Indeed a comparison of 
figures [6l and 1 1 1 1 shows that a cluster can have a smaller value 
of rhm than another cluster but still have a larger value of 
rhi- In figure [T4l I show the ratio of r^i to rhm for all of my 
cluster models. The ratio varies considerably both with time 
and with the initial conditions used. The ratio does seem 
to attain a fairly stable value for each individual cluster, 
indicating that Vhi is not completely independent of rhm- 
Since this value always lies somewhere between 0.35 and 0.7 
it is, in principle, possible to use Vhi to estimate the value of 
rhm to within a factor of ~ 2. It is not possible, however, to 
use Rhi,br to estimate the relative size of two clusters unless 
(at minimum) their metallicities, initial concentrations, and 
ages are known. 



6 CONCLUSIONS 

I have analyzed a suite of Monte Carlo globular cluster sim- 
ulations to determine the origin of the difference in half-light 
radii between blue and red GCs. I find that, provided they 
are mass-segregated, the value of rhi is larger in blue than 
in red GCs simply due to differences in the luminosity func- 
tion between metal-poor and metal rich stellar populations. 
Depending on the initial conditions this effect is sufficient 
to explain the observed differences between rhi in blue and 
red GCs. I find that it is also possible for blue clusters to be 
physically larger than red clusters as measured by their half- 
mass radii. However this only occurs if there are significant 
differences in the number and properties of BHs between 
metal-poor and metal-rich stellar populations. A difference 
in rhm can certainly enhance the difference in rhi but it 
is not a necessary condition for such a difference to exist. 
Furthermore, the enhancement in rhi in blue GCs when a 
there is significant difference in the BH population is larger 
than observed. This leads me to conclude that blue and red 
GCs probably do not have significant differences in their BH 
populations but further simulations and more careful com- 
parisons with observations will be necessary to confirm this. 
Further simulations are also necessary to determine if it is 
number or mass that determines whether BHs can cause a 
GC to expand. Finally, I find that Rhi,br does not directly 
predict Rhm,br so differences in rhi cannot be used to infer 
differences in size between clusters of different metallicities. 
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